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We determine under what conditions Scalar Tensor cosmologies predict an expansion rate which 
is reduced as compared to the standard General Relativity case. We show that ST theories with 
a single matter sector typically predict an enchanced Hubble rate in the past, as a consequence of 
the requirement of an attractive fixed point towards General Relativity at late times. Instead, when 
additional matter sectors with different conformal factors are added, the late time convergence to 
General Relativity is mantained and at the same time a reduced expansion rate in the past can be 
driven. For suitable choices of the parameters which govern the scalar field evolution, a sizeable 
reduction (up to about 2 orders of magnitude) of the Hubble rate prior to Big Bang Nucleosynthesis 
can be obtained. We then discuss the impact of these cosmological models on the relic abundance 
of dark matter is minimal Supergravity models: we show that the cosmologically allowed regions in 
parameter space are significantly enlarged, implying a change in the potential reach of LHC on the 
neutralino phenomenology. 

PACS numbers: 95.35.+d,95.36.+x,98.80.-k,04.50.+h,98.80.Cq 



I. CHANGING THE EXPANSION RATE IN 
THE PAST 

In a standard flat FRW universe described by GR, the 
expansion rate of the universe, Hgr = a/ a, is set by the 
total energy density, /5 to t, according to the Friedmann 
law, 

where M p is the Planck mass, related to the Newton con- 
stant by M p = (87rG) _1 / 2 . If the total energy density 
is dominated by relativistic degrees of freedom, the ex- 
pansion rate is related to the temperature through the 
relation 

H GR ~ 1.66 gl /2 ^~, (2) 

with g» the effective number of relativistic degrees of free- 
dom (see for instance [1]). 
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In order to modify the above H-T relation, one can do 
one (or more) of the following: 

1) change the number of relativistic d.o.f.'s, g*; 

2) consider a p to t not dominated by relativistic d.o.f.'s; 

3) consider a modification of GR in which an effective 
Planck mass, different from M p appears in (2). 

One example of a scenario of the first type is obtained 
by adding N extra light neutrino families to the standard 
model, which increases by 7/4 TV. 

The second situation is realized e.g. in the so called 
"kination" scenario [2], where the energy density at a 
certain epoch is dominated by the kinetic energy of a 
scalar field. Since the kinetic energy redshifts as pkin ~ 
a~ 6 , it will eventually become subdominant with respect 
to radiation (p ra d ~ a -4 ). As long as pkin dominates, the 
expansion rate is bigger than in the "standard" scenario 
where there is no scalar field and the same amount of 
radiation. 

In this paper we will consider scenarios of the third 
type, where the expansion rate is modified by changing 
the effective gravitational coupling. This can be realized 
in a fully covariant way in ST theories [3]. We will con- 
sider the class of ST theories which can be defined by the 
following action [4], 

S = S g + J2S h (3) 
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where S g is the gravitational part, given by the sum of 
the Einstcin-Hilbcrt and the scalar field actions, 



Ml 



Sg = I d^xj—g R + g^d^d u ip - j^V(p) 

'(4) 

where V(ip) can be either a true potential or a (Einstein 
frame) cosmological constant, V((p) = Vq- The Si's are 
the actions for separate "matter" sectors 



(5) 



with indicating a generic field of the i-th matter sec- 
tor, coupled to the metric A 2 {p>)g^ v . The actions Si 
are constructed starting from the Minkowski actions of 
Quantum Field Theory, for instance the SM or the MSSM 
ones, by substituting the flat metric everywhere with 
A 2 (<p)g^. 

The emergence of such a structure, with different con- 
formal factors A 2 for the various sectors can be motivated 
in extra-dimensional models, assuming that the two sec- 
tors live in different portions of the extra-dimensional 
space. A similar structure, leading to more dark mat- 
ter species, each with a different conformal factor, was 
considered in [5]. 

We consider a flat FRW space-time 

ds 2 = dt 2 - a 2 {t) dl 2 , 

where the matter energy-momentum tensors, = 

2(—g)~ 1 / 2 SSi/8g ,J ' 1 ' admit the perfect-fluid representa- 
tion 



T U = (fit +Pi) u V u v ~Pi 9fiv , 

with g^ v u^u u = 1. 

The cosmological equations then take the form 
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where the coupling functions ct^ are given by 

d\og Ai 
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dp 
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The Bianchi identity holds for each matter sector sepa- 
rately, and reads, 



d(pi a 3 ) +p t da 3 = (pi - 3 pi) a 3 d\ogAi(<p), 
implying that the energy densities scale as 



M<f) 
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(11) 



(12) 



with Wi = Pi/ Pi the equation of state associated to the 
i-th energy density (assuming Wi is constant). 



A. GR as a fixed point 

To start, consider the case of a single matter sector, 
Sm- In order to compare the ST case with the GR 
one of Eqs. (1, 2), it is convenient to Weyl— transform 
to the so-called Jordan Frame (JF), where the energy- 
momentum tensor is covariantly conserved. The trans- 
formation amounts to a rescaling of the metric according 
to 



g^iv = A 2 M (ip)g^ , 



(13) 



keeping the comoving spatial coordinates and the con- 
formal time drj = dt/a fixed [6], The JF matter energy- 
momentum tensor, T™ = 2(—g)~ 1 ^ 2 8SM /Sg^ v , is re- 
lated to that in eq. (6) by TM, = Aj 2 T^l, so that energy 
density and pressure transform as 
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while the cosmic time transforms as dt = Amdt. One can 
easily verify that the above defined quantities satisfy the 
usual Bianchi identity, that is Eq. (11) with vanishing 
RHS, and that, as a consequence, pu ~ a _3 ( 1 + u,M ). The 
expansion rate, H$t = dlogd/dt, is given by 



oc M {<p)p' a 



A M (p>) 



(15) 



where we have defined um according to Eq. (10), and 
(•)' = d(-)/dloga. Using (15) and (14) in (8), we obtain 
the Fricdmann equation in the ST theory, 
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where V = A M V. Comparing to Eq.(l), we see that 
apart from the extra contribution to p to t from the scalar 
field potential, the ST Friedmann equation differs from 
the standard one of GR by the presence of an effective, 
field-dependent Planck mass, 



A 2 M (P) (1 
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where the last equality holds with very good approxima- 
tions for all the choices of Ai functions considered in the 
present paper. 

If the conformal factor A 2 M (p) is constant, then the full 
action S 3 + Sm is just that of GR (with M p = M*/Am) 
plus a minimally coupled scalar field. Therefore, the cou- 
pling function ajj, defined according to Eq. (10), mea- 
sures the "distance" from GR of the ST theory, a a/ = 
being the GR limit. Changing Am, and, therefore, 
changing the effective Planck mass, opens the way to 
a modification of the standard relation between H and 
p, or T. In order to study the evolution of A M (p>), one 
should come back to Eq. (9). Considering an initial epoch 
deeply inside radiation domination, we can neglect the 
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contribution from the potential on the RHS. The other 
contribution, the trace of the energy-momentum tensor 
(Pm — 3 pm) is zero for fully rclativistic components but 
turns on to positive values each time the temperature 
drops below the mass threshold of a particle in the ther- 
mal bath. Assuming a mass spectrum - e.g that of the 
SM, or of the MSSM - one finds that this effect is effec- 
tive enough to drive the scalar field evolution even in the 
radiation domination era [7]. 

The key point to notice is that if there is a field value, 
tpo, such that cxm{<Po) = 0, this is a fixed point of the field 
evolution [8, 9]. Moreover, if a' M is positive (negative) 
the fixed point is attractive (repulsive). Since olm = 
corresponds to the GR limit, we see that GR is a - 
possibly attractive - fixed point configuration. 

The impact on the DM relic abundance of a scenario 
based on this mechanism of attraction towards GR was 
considered in [7, 10]. Regardless of the particular form of 
the j4m(^) function, the requirement that an attractive 
fixed point towards GR exists implies that the effective 
Planck mass in the past was not smaller than today, that 
is to say that, at a certain temperature T, for instance 
at DM freeze out, the universe was expanding not more 
slowly than in the standard GR case. This is easy to 
understand, since the past values of tp, and then of Am(<p) 
are all such that 



log 



A M (tp) 
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dx olm{x) > , 
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with tp between tpo and the next fixed point. Therefore, 
according to Eq. (17), the ratio between Hst and Hqr, 



Hst 

H-GR 



(19) 



can only decrease in time. Another way of seeing this, is 
by noticing that the RHS of the field equation (9), is pro- 
portional to the field derivative of the effective potential 
Vtff = Pm + y , where the field dependence of pm is given 
by Eq. (12). Then, neglecting again V, the field evolu- 
tion will tend towards minimizing Am(v) (if %< 1/3), 
therefore minimizing Hst/ Hqr. 

Eq. (19) is obtained under the same approximation 
used in Eq. (17), that is by neglecting the scalar field 
contribution to the total energy density, and assuming 
the same matter content in the ST and GR cases. Notice 
that these approximations are far from mandatory, and 
we only use them here in order to illustrate how the fixed 
point mechanism works. The numerical analysis we will 
present in the following were indeed obtained using the 
full expressions, such as Eq. (16). Finally, in order to 
identify the fixed point with GR, we have to impose 
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B. Lowering H in the past 

All the mechanisms discussed so far (i.e. adding rela- 
tivistic d.o.f.'s, the kination scenario, or ST theories with 
a single matter sector) give a faster expansion of the uni- 
verse in the past w.r.t. the standard case. In the case 
of ST theories with a single matter sector, we have just 
seen that this comes as a consequence of the requirement 
of an attractive fixed point towards GR. In this subsec- 
tion, we will show that adding more matter sectors, with 
different conformal factors, allows us to keep the desir- 
able property of late time convergence to GR and, at the 
same time, to have a lower expansion rate in the past. 
To illustrate this point, we will consider just two matter 
sectors, a "visible" one, containing the SM or one of its 
extensions, and a "hidden" one. The full action is then 
given by 



S — Sn ~\~ S v ~\~ Sfi . 



(21) 



where the two matter actions S v and Sh have two differ- 
ent conformal functions A v (tp) and Ah(tp). The discus- 
sion follows quite closely that of the previous subsection. 
The only subtle point is to notice that, if A v [tp) ^ Ah(tp) 
there is no Weyl transformation that gives covariantly 
conserved energy-momentum tensors both for the visible 
and for the hidden sector. Since particle masses, reaction 
rates and so on, are computed in terms of parameters of 
the "visible" action, the transformation to perform in or- 
der to compare with the standard GR case is the one 
leading to a conserved TJJ„, that is [4] 



g» v = A 2 v (tp)g^ 



(22) 



and so on. As a consequence, the expansion rate in this 
case is given by 



H ST — 
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where 
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In order to study the existence of a fixed point, it is still 
convenient to revert to the Einstein Frame field equation, 
Eq. (9). The RHS, is now given by the field derivative of 
the effective potential 



V eff = Pv + Ph + V , 



(24) 



with the field-dependence of p Vi h given by Eq. (12). The 
condition to have a fixed point is then 



on (1 - 3Wi) ft + V = , 



(25) 



i—v.h 
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while, asking that the fixed point is stable implies 

(a'i (1 - 3w i )p i + a\ (1 - Zw l f Pl + V") > . 



i—v,h 



(26) 

From Eq. (23) we see that, away from the fixed point, 
Hst is lower than the one obtained in GR with the same 
matter content but frozen scalar fields if 



d 2 ( A 2 ( ^ l + Ph/Pv\sT 



<0, (27) 



where, again, we have assumed that the second fraction 
in Eq. (23) is approximately one, and we have neglected 
the contribution from the scalar potential. 

As an example, we now consider Ai functions of the 
form 



(28) 



Neglecting again the potential, we see that the fixed point 
condition, Eq. (25), is solved by the symmetric point (p = 
0. The stability condition, Eq. (26), translates into 



J2 bi(l-3wi) Pi >0, 



(29) 



i—v.h 



which, according to Eq. (27), is compatible with a lower 
Hst / Hqr outside the fixed point {i.e. in the past), if 



b v < 0, 



(30) 



where we have assumed ph <C p v close to the fixed point, 
since we are interested in a physical situation in which 
most of the dark matter lives in the "visible" sector (as 
in the MSSM). 



C. Numerical examples 

To be implemented in a sensible cosmological model, 
the previously discussed mechanism for lowering the ex- 
pansion rate in the past has to respect the severe bound 
coming from BBN, namely [11] 



\Hst — Hgr\ 
Hst 



< 10% 



at BBN . 



(31) 



We now show in a few examples that indeed the 
bound (31) can be satisfied even by points of the param- 
eters space {b v ,bh,<pin) giving rise to a pre-BBN value 
of the ratio (19) as low as 10 -3 . Such important devia- 
tions from standard cosmology are allowed in the present 
scenario by the effectiveness of the GR fixed point. 

In order to numerically solve Eqs. (7-9) we need the 
equation of state parameters 
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i = v,h (32) 



where y = logo and y 1 refers to the equivalence in the 
visible (i = v) and hidden (i = h) sectors respectively. 
The functions Ii(y) are given by [7] 
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where yp i = — logmpjTo, mp i and gp i are the masses 
and the relativistic degrees of freedom of the particles 
and T = 2.73 K ~ 2.35 x 10" 13 GcV is the current 
temperature of the Universe. 

As mass thresholds for the visible sector, namely yp v , 
we use the one given by a MSSM-like mass spectrum 
[36] while the equivalence time y° eq has been computed 
according to [12]. The analogous quantities for the hid- 
den sector are free parameters of the theory that nev- 
ertheless do not have any drastic impact on the final 
result. The only significant assumption we do is that 
lim a ^(l — 3u>h(y)) — 1 for y <C yBBN ■ In such a way 
the contribution of the hidden sector to the RHS of the 
scalar field equation before BBN can be dominant w.r.t. 
the one of the visible sector. The equation of state pa- 
rameter Wh is plotted in Fig. (1) for three different values 
of y^ q . In Fig. (2) instead, for a given value of y^ q we plot 
Wh for three different choices of mass thresholds in the 
hidden sector. 

We can now integrate the equation of motion for the 
scalar field. Fixing as the initial condition ifi n = 1 (in 
Planck units), we show in Fig. (3) the behavior of tp as 
a function of y for three different choices of parameters 
b v and bh- For the same choices we plot in Figs. (4) 
and (5) the evolution of the ratio Hg T j Hq R and of the 
function a v respectively As anticipated, an agreement 
with the BBN bound is achieved even by solutions with 
a pre-BBN value of the ratio H st /Hq R of order 10~ 3 . 

The scalar field dynamics can be qualitatively under- 
stood by having in mind the following rough estimation 
for its first derivative 



if oc 



2jaj(l - 3w t )p l 



(34) 



where each ai is weighted by the the corresponding func- 
tion (1 — 3vii)pi. If both ai are positive, the scalar field 
is driven toward the fixed point 4> — 0. Analogously, 
negative values of the couplings ai lead to a run-away 
behavior of ip. The present scenario corresponds to the 
case in which (before BBN) a v < and ah > 0. There- 
fore, the interplay between the contributions of the hid- 
den and visible sectors to the RHS of Eq. (34) becomes 
relevant. If at early time (1 — 3wh)ph always dominates, 
then the only effect of having a a v < is to realize the ini- 
tial condition Hg T / Hq R < 1 (see Subsection IB). This 
is the case in Models 1 and 2 in Fig. (4). However, if 
at early time there is a short period where the visible 



5 




FIG. 1: Four different combinations of the quantity (1 — 3u>i) as a function of y. For all of them a MSSM-like mass spectrum 
is assumed (heaviest mass 1 TeV). The solid [green] line corresponds to the visible sector; the dot-dashed [red], dashed [blue] 
and dotted [black] lines represent a hidden sector with a hidden matter-hidden radiation equivalence at 10 GeV, 1 GeV and 
0.1 GeV respectively. 




FIG. 2: Four different combinations of the quantity (1 — 3wi) as a function of y. The solid [green] line corresponds to the 
visible sector with a MSSM-like mass spectrum (heaviest mass 1 TeV). The dot-dashed [red], dashed [blue] and dotted [black] 
lines represent a hidden sector with three different mass spectra: heaviest mass 0.1 TeV, 1 TeV and 10 TeV respectively. In 
the hidden sector the equivalence temperature is fixed at 1 GeV. 



sector contribution (1 — 3w v )p v dominates, then the ra- 
tio Hg T / Hq R decreases from its initial value until when 
(1 — 3w v ) ~ 1 and, as a consequence, the hidden sector 
contribution becomes dominant. This happens in Model 
3 of Fig. (4) where the parameter b v has been tuned close 
to -1/2 [37]. 



choosing A = b v ~ — 1/2 we find 

A b „ ~ 2|(1 + b v )b v \ 10 3 - 5 x 10 2 . (35) 

Therefore, configurations with a small initial value of 
Hg T / Hq R are very sensitive to the parameter b v . 



Let us conclude this subsection with an estimation of 
the dependence of our results from the parameters b v . 
According to [13], the level of fine-tuning A\ on a pa- 
rameter A needed to get the required value of an observ- 
able O is given by A A = \{\/0){dO/d\)\ . Requiring 
at early time (when ip ~ 1) 0= H"^ t /Hq R ~ 10 -3 and 



II. IMPLICATIONS FOR DARK MATTER IN 
THE CMSSM 

A modification of the Hubble rate at early times has 
impact on the formation of dark matter as a thermal 
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FIG. 3: Evolution of the scalar field (p with y for three different choices of the parameters (b v ,bh)- The dark solid [red], dashed 
[blue] and light solid [green] lines correspond respectively to (-0.2,5) [Model 1], (-0.4, 15) [Model 2] and (-0.521,50) [Model 
3]. BBN occurs at y ~ -22. 
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FIG. 4: Ratio of the ST and GR Hubble rates squared Hg T / Hq R as a function of y for three different choices of parameters 
(b v ,bh)- The dark solid [red], dashed [blue] and light solid [green] lines correspond respectively to (—0.2, 5)[Model 1] , (—0.4, 15) 
[Model 2] and (-0.521,50) [Model 3]. BBN occurs at y ~ -22. 



relic, if the particle freeze-out occurs during the period 
of modification of the expansion rate. ST cosmologies 
with a Hubble rate increased with respect to the GR 
case have been discussed in Rcfs. [7, 10, 14], where the 
effect on the decoupling of a cold relic was discussed 
and bounds on the amount of increase of the Hubble 
rate prior to Big Bang Nucleosynthesis have been de- 
rived from the indirect detection signals of dark matter in 
our Galaxy. For cosmological models with an enhanced 
Hubble rate, the decoupling is anticipated, and the re- 
quired amount of cold dark matter is obtained for larger 
annihilation cross-sections: this, in turn, translates into 
larger indirect detection rates, which depend directly on 
the annihilation process. In Refs. [10, 14] we discussed 



how low-energy antiprotons and gamma-rays fluxes from 
the galactic center can pose limits on the admissible en- 
chancement of the pre-BBN Hubble rate. We showed 
that these limits may be severe: for dark matter parti- 
cles lighter than about a few hundred GeV antiprotons 
set the most important limits, which arc quite strong for 
dark matter masses below 100 GeV. For heavier parti- 
cles, gamma-rays are more instrumental in determining 
significant bounds. Further recent considerations on the 
effect of cosmologies with modified Hubble rate are dis- 
cussed in Refs. [15, 16, 17, 18, 19, 20]. 

In the case of the cosmological models which predict a 
reduced Hubble rate, the situation is opposite: a smaller 
expansion rate implies that the cold relic particle remains 
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FIG. 5: Evolution of the coupling function a v with y for three different choices of parameters {b v ,bh)- The dark solid [red], 
dashed [blue] and light solid [green] lines correspond respectively to (—0.2,5) [Model 1], (—0.4, 15) [Model 2] and (—0.521,50) 
[Model 3]. The BBN occurs at y ~ -22 



in equilibrium for a longer time in the early Universe, 
and, as a consequence, its relic abundance turns out to 
be smaller than the one obtained in GR. In this case, the 
required amount of dark matter is obtained for smaller 
annihilation cross sections, and therefore indirect detec- 
tion signals are depressed as compared to the standard 
GR case: as a consequence, no relevant bounds on the 
pre-BBN expansion rate can be set. On the other hand, 
for those particle physics models which typically predict 
large values for the relic abundance of the dark matter 
candidate, this class of ST cosmologies may have an im- 
portant impact in the selection of the regions in param- 
eter space which arc cosmologically allowed. 

A typical and noticeable case where the relic abun- 
dance constraint is very strong is offered by minimal 
SUGRA models, where the neutralino is the dark matter 
candidate and its relic abundance easily turns out to be 
very large, in excess of the cosmological bound provided 
by WMAP [12]: 

0.092 < Ocdm^ 2 < 0.124 (36) 

Large sectors of the supersymmetric parameter space are 
excluded by this bound. A reduction of the expansion 
rate will therefore have a crucial impact on the allowed 
regions in parameter space, which are therefore enlarged. 
The potential reach of accelerators like the Large Hadron 
Collider (LHC) or the International Linear Collider (ILC) 
on the search of supersymmetry may therefore be affected 
by this broadening of the allowed parameter space, espe- 
cially for the interesting situation of looking for super- 
symmetric configurations able to fully explain the dark 
matter problem. 

We have therefore studied how the allowed parame- 
ter space of minimal SUGRA changes in ST cosmologies 
with a reduced Hubble rate. We have used a cosmo- 
logical model of the type of Model 3 discussed in the 



previous Section and depicted in Fig. 4. For the calcu- 
lations of the neutralino relic density we have used the 
DarkSUSYpackage [21], with an interface to ISAJET 7.69 
[22] for the minimal SUGRA parameter space determi- 
nation, with two major modifications. First, the relic 
density is obtained by the implementation of a numeri- 
cal solution of a modified Boltzmann equation which in- 
cludes the reduced Hubble rate evolution (similar to the 
method used in Refs. [7, 10, 14] for the enhanced case). 
Second, the NNLO contributions to the Standard Model 
branching ratio of the BR[£? — > X s j] have been recently 
determined [23]: the updated result, which we use here, 
is BRL9 -> X s -f] SM = (3.15 ± 0.23) x lO" 4 (£ 7 > 1.6 
GeV). In order to implement the NNLO SM result with 
the supersymmetric contribution, which are known up 
to the NLO [24] , we have used the following approximate 
expression, which is suitable when the beyond-standard- 
model (BSM) corrections are small [23, 25, 26]: 

[BRjtheory X 10 4 = 3.15 (37) 

-1.9x(w[cf] + ^)5 BSM [C«]) 

where CfVo) are LO (j = 0) and NLO (j = 1) 
Wilson coefficients. For the matching scale /io (which 
should be taken as /j,q = 2Mw ~ 160 GeV) we use 
/io = TOt(mt, pole) = 163.7 GeV, and we use a top-quark 
pole mass m t , P oi c = 171.4 GeV [27] and a s (M z ) = 0.1189 
[28]. The theoretical calculation in Eq. (37) is compared 
to the current world average of the experimental deter- 
mination [29]: 

[ BR L P x 104 = f 3 - 55 ± °- 24 in'?n ± °- 03 ) ( 38 ) 
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FIG. 6: Regions in the (m^^o) parameter space where the 
neutralino relic abundance falls in the cosmological interval 
for cold dark matter obtained by WMAP, for tan f3 = 10, 
Aq — and positive fj,. In the bulk region, the lower [black] 
points refer to GR cosmology, while the upper [red] points 
stand for a ST cosmology with a reduced Hubble rate. The 
shaded areas are forbidden by theoretical arguments and ex- 
perimental bounds. The two curves are indicative of the reach 
for 100 fb^ 1 of the LHC [33, 34] and of the ILC at y/s = 1000 
GeV energy [33, 35]. 

The estimated error of the theoretical calculation of the 
Standard Model contribution is ±0.23 x 10" 4 [23, 25]. For 
the theoretical beyond-SM correction we have assumed 
an error of the same size. Adding all experimental and 
theoretical errors we get the following 2er interval for the 
branching ratio: 

2.71 x 1CT 4 < BR[B -> X sl ] < 4.39 x 1(T 4 (39) 

Supcrsymmetric models for which the theoretical calcu- 
lation in Eq. (37) is outside this interval are considered 
in disagreement with the experimental result. 

A. Low tan/3 

As a first example, we scan the universal gaugino mass 
m-y/2 and soft scalar-mass too parameters of minimal 
SUGRA for a low value of the tan/3 parameter (tan/3 
is defined, as usual, as the ratio of the two Higgs vacuum 
expectation values V2 and v\, where vi (vi) gives mass 
to the top(down)-type fermions) and a vanishing univer- 
sal trilincar coupling Aq. The higgs-mixing parameter [i, 
derived by rcnormalization group equation (RGE) evolu- 
tion and electro-weak symmetry breaking (EWSB) con- 
ditions, is taken as positive. Our choice of parameters is 



here: 

tan/3 = 10 sgn(^) = + A = (40) 

Fig. 6 shows the result in the plane (toi/2, too), both for 
the standard GR case and for the ST case of Model 3. 
Shaded areas denote regions which are excluded cither 
by theoretical arguments or by experimental constraints 
on higgs and supcrsymmetry searches, as well as supcr- 
symmetric contributions to rare processes, namely to the 
BR[£> — > X^] and to the muon anomalous magnetic mo- 
ment (g — 2)^. More specifically, the upper wedge refers 
to the non-occurrence of the radiative EWSB and the 
lower-right area to the occurrence of a stau LSP. The 
low-mi/2 vertical band is excluded by the quoted exper- 
imental bounds. 

The sector of the supcrsymmetric parameter space 
which provides LSP neutralinos with a relic abundance 
in the cosmological range of Eq. (36) are denoted by the 
open circles: in the so-called "bulk region" (low values 
of both TO1/2 and too), the lower [black] points fulfill the 
density constraint in the standard GR cosmology, while 
the upper [red] points in the modified ST cosmology with 
reduced Hubble rate. In the region above the points, 
the neutralino relic abundance exceeds the cosmological 
bounds, and therefore refers to supcrsymmetric configu- 
rations which are excluded by cosmology. Fig. 6 shows 
that in our modified cosmological scenario, the allowed 
regions in parameter space are enlarged (the relic density 
has been decreased 1.4 times to 4.4 times compared to the 
standard case) and those which refer to dominant neu- 
tralino dark matter are shifted towards larger values of 
the supcrsymmetric parameters toi/2 and Too- The bulk 
region now occurs for values of mo larger by a factor of 2 
(while the bulk region for the GR case now refers to cos- 
mologically subdominant neutralinos). Nevertheless, this 
sector of the parameter space is already mostly excluded 
by accelerator searches. In the coannihilation channel, 
which extends for low values of the ratio Too/to^, along 
the boundary of the stau excluded region, the change is 
more dramatic: this coannihilation region, which appears 
to be fully exploitable at the LHC, now extends towards 
larger values of mi/2, beyond the estimated LHC reach. 

In the cosmologically allowed region of large too , where 
a gaugino-mixing becomes possible and therefore the 
neutralino can efficiently annihilate and provide an ac- 
ceptable relic abundance (a mechanism discussed in Ref. 
[30] and lately dubbed as "focus point region" in Ref. 
[31]), the effect of reducing the Hubble rate translates 
into a slight lowering of the cosmologically relevant re- 
gion, with no drastic phenomenological effect in this case. 

B. Large tan /3 

As a second example, we consider the case of large 
values of tan/3. We show two cases, one which refers to 
a negative /i parameter: 

tan/3 = 45 sgn(/i) = - A a = (41) 
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FIG. 7: Regions in the (toi/2,T71o) parameter space where the neutralino relic abundance falls in the cosmological interval for 
cold dark matter obtained by WMAP. The left panel refers to tan/3 = 45, Ao — and negative fj,. The right panel is obtained 
for tan (3 = 53, Aq = and positive fx. Notations are as in Fig. 6. 



and one to a positive value of /Lt: 

tan f3 = 53 sgn(» = + A = (42) 

The results are shown in Fig. 7. In this case the change 
in the cosmological scenario is more relevant, not only in 
the coannihilation channel, but also in the "funnel" region 
[ ] which occurs for intermediate values of the ratio 
m /m 1 / 2 - In the GR case, almost the full cosmologi- 
cally allowed parameter space may be explored by the 
LHC. When the ST cosmology is considered, the funnel 
region dramatically extends towards large values of tuq 
and mi/2 an d goes well beyond the accelerators reach. 
Also the coannihilation region now extends to values of 
7?7i/2 well in excess of 2 TeV. In the case of the posi- 
tive fi reported in the right panel of Fig. 7, also the 
focus-point region shows a deviation from the GR case, 
and is shifted towards lower values of the mo parame- 
ter. In summary, for these large values of tan j3 the reach 
of LHC on the cosmologically relevant configurations is 
less strong than in the case of GR: a discover of super- 
symmetry will be more likely related to a subdominant 
neutralino dark matter. 



III. CONCLUSIONS 

We have discussed Scalar Tensor cosmologies by deter- 
mining under what conditions these theories can predict 
an expansion rate which is reduced as compared to the 
standard General Relativity case. We showed that in the 



case of ST theories with a single matter sector, the typ- 
ical behaviour is an enchancement of the Hubble rate in 
the past: this arises as a consequence of the requirement 
of an attractive fixed point towards GR at late times. 
Instead, when additional matter sectors, with different 
conformal factors, are added, we can mantain the desir- 
able property of late time convergence to GR and, at the 
same time, obtain a reduced expansion rate in the past. 
We showed that, for suitable choices of the parameters 
which govern the scalar field evolution, a sizeable reduc- 
tion (up to about 2 orders of magnitude) of the Hubble 
rate prior to Big Bag Nucleosynthesis can be obtained. 
Large reductions come along with some fine-tuning on 
the scalar field parameters, while a milder decrease oc- 
curs without tuning problems. 

We have then applied the results obtained on the re- 
duction of the early-time Hubble rate to the formation 
of dark matter and the determination of its relic abun- 
dance. If the dark matter decouples during the period of 
Hubble-rate reduction, the relic abundance turns out to 
be reduced as compared to the standard GR case. This 
has therefore impact on the determination of the cos- 
mologically allowed parameter space in minimal SUGRA 
models, where, in large portions of the parameter space 
and for the GR case, the neutralino relic abundance is 
large and in excess of the WMAP bound. We have there- 
fore explicitely shown what are the modifications to the 
minimal SUGRA allowed parameter space when ST cos- 
mologies with a reduced Hubble rate are considered and 
we have quantified the effect in view of the reach of LHC 
and ILC on the searches for supersymmetry at future ac- 
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cclcrators. These modifications move the cosmologically 
relevant regions up to a few TeV for the mwa param- 
eters, since they significantly extend the coannihilation 
corridor and the funnel region which occurs at large val- 
ues of tan/3. 
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